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Abstract 

We compute the three-dimensional effective action for the minimal supersymmetric stan- 
dard model, which describes the light modes of the theory near the finite-temperature elec- 
troweak phase transition, keeping the one-loop corrections from the third generation quarks 
and squarks. Using the lattice results of Kajantie et al. for the phase transition in the same 
class of 3-D models, we find that the strength of the phase transition is sufficient for elec- 
troweak baryogenesis, in much broader regions of parameter space than have been indicated 
by purely perturbative analyses. In particular we find that, while small values of tan/3 are 
favored, positive results persist even for arbitrarily large values of tan/3 if the mass of the 
A boson is between 40 and 120 GeV, a region of parameters which has not been previously 
identified as being favorable for electroweak baryogenesis. 



1 Introduction 



One of the fundamental questions in nature is the origin of the asymmetry of matter over 
antimatter in the universe. Although explanations abound, one of the most interesting pos- 
sibilities is that the baryon asymmetry was created during the electroweak phase transition 
(EWPT) [Q, using new physics at sufficiently low energies to be verifiable in anticipated 
experiments like LEP-II or the Large Hadron Collider. Although the EWPT is a first order 
transition in the standard model, it is too weakly so to fulfill Sakharov's out-of-thermal- 
equilibrium requirement for generating baryons: any asymmetry created during the EWPT 
would be quickly erased afterwards by residual sphaleron interactions in the broken phase of 
the SU(2)l x U(l) gauge theory |fj]. Moreover it appears that the standard model has too 
little CP violation for electroweak baryogenesis || . 

It is therefore interesting to find out whether a more strongly first order EWPT is possible 
in extensions of the Standard Model, a prime example being its minimal supersymmetric 
extension, the MSSM, shown to be suitable for electroweak baryogenesis in ref. [|J. The 
phase transition has been studied in this model by means of the one-loop finite-temperature 
effective potential with the result that there exist some regions of parameter space 

where electroweak baryogenesis is possible. However the perturbative approach should be 
viewed with skepticism because at finite temperature it becomes infrared divergent for small 
values of the Higgs field, which can be crucial for determining the critical temperature and 
Higgs field VEV at the phase transition. 

To deal with the breakdown of perturbation theory, Kajantie et al. |2j have studied the 
EWPT of the Standard Model on the lattice. As they have emphasized however, almost any 
extension of the Standard Model can be reduced to an effective three-dimensional theory of 
one Higgs doublet interacting with SU(2) gauge bosons, by integrating out all the modes 
with thermal masses larger than those of the longitudinal gauge bosons ||. The beauty of 
this approach is that the effective theory need only be numerically studied once; after that it 
is simply a matter of matching the parameters of the fundamental theory onto this effective 
theory, which can be reliably done using perturbation theory. 

In this paper we construct the effective 3-D Lagrangian corresponding to the MSSM 
at finite temperature, keeping the dominant effects proportional to the top quark Yukawa 
coupling. At the critical temperature it has the simple form 

£ 3 = \(di - \9 ? >t ■ W)<$>\ 2 + m 2 &§ + A 3 ($ f $ ) 2 + \F l3 F ir (1) 



1 



where the couplings A3 and #3 depend, ultimately, on physical parameters of the MSSM such 
as tan /3, Higgs boson masses, and squark masses. The criterion from lattice studies || for 
preserving any baryon asymmetry created during the EWPT is 

^ < 0.04. (2) 

#3 

We find that this bound is satisfied in a significant fraction of the MSSM parameter space, 
in contrast with previous studies based on a purely perturbative approach || . Although the 
most recent perturbative investigations obtained more positive results by considering nega- 
tive values of the squark mass parameter mfj |7J or small values of the right-handed squark 
||, in the present work we find that such choices, while compatible with electroweak baryo- 
genesis, are not particularly favored and represent only a small fraction of the total volume 
of baryogenesis-allowed parameter space. The quantities to which our results turn out to 
be most sensitive are the ratio of Higgs VEV's, (H 2 ) / (H\) = tan/3, the mass of the pseu- 
doscalar Higgs boson A and the soft supersymmetry breaking squark mixing parameters. 
We will show that the allowed regions can be characterized roughly by 0.5 < tan/3 < 2 
and m^o unrestricted, or m^o between 40 and 120 GeV for arbitrarily large tan/3, and no 
special restrictions on the other parameters except that the largest portion of the allowed 
space corresponds to large squark mixing parameters. This would appear to be a much 
less constrained situation than was previously believed to exist for the MSSM as regards 
baryogenesis. 

The most important interactions affecting the strength of the phase transition are those 
involving the largest couplings to the Higgs field, namely the top quark Yukawa coupling y t . 
It is conceivable that tan/3 ~ m t /mb, in which case the bottom quark Yukawa coupling y b 
would also be large. Then the relevant part of the MSSM Lagrangian, including the neutral 
sector of the two Higgs fields (so the below are not doublets) and the third generation 
quarks and squarks, can be written in Euclidean space as 

£ tree = £ (| Ml | 2 + m J 2 |^| 2 )+m^(^2 + h.c.) + ^4^(l^| 2 -|^| 2 ) 2 

i=l,2 8 

+ yti L H 2 t R + h.c. + y 2 \H 2 \ 2 (\t L \ 2 + |£r| 2 ) + yJXfiHi + A t H 2 )t R + h.c. 
+ ythHibR + h.c. +yl\H 1 \ 2 (\b L \ 2 + \b R \ 2 )+y b b* L {fiH 2 + A b H 1 )b R + h.c. 

+ (b 2 (m 2 - \h\ 2 ) - y 2 (m 2 + \b L \ 2 ) + y 2 (2\t R \ 2 - \i R \ 2 )) {m 2 - \h 2 \ 2 ) 

+ ™ 2 Q (\i L \ 2 + \b L \ 2 ) +m 2 u \t R \ 2 + m 2 D \b R \ 2 . (3) 
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Figure 1: The 1PI diagrams needed for the scalar 2-point function. Thin line in the loop 
represents quarks and the heavy lines squarks. Indices labeling external legs refer to doublet; 
when index is a letter either 1 or 2 is allowed. The example shown corresponds to top quark 
or squark in the loop; for bottom sector reverse 1 and 2. 

where the mf's for % — 1,2,3 contain both the soft-breaking and the supersymmetric con- 
tributions. In the present work we will consider only moderatly large values for tan/3, so 
that yb <C y t in the numerical analysis; the terms of order yb are nevertheless displayed, to 
allow for future investigation of the large tan/3 regime. However even if yj, <C y t , the bottom 
squarks can still be relevant because of certain loop diagrams proportional to masses, for 
which they contribute competetively with the top squarks if they are sufficiently heavy, and 
these effects we do include throughout. 

To apply the lattice gauge theory bound (|2]) one must carry out two steps ||. First, 
integrate out all the heavy degrees of freedom in the finite-temperature theory at the phase 
transition. This means everything except for a single light linear combination of the Higgs 
fields, and the transverse gauge bosons, resulting in the three-dimensional effective action 
(H) for these fields. Second, renormalize the same theory at zero temperature so as to express 
the parameters appearing in £3 as functions of physical observables, such as particle masses. 
In both steps we compute the corrections due to third generation quarks and squarks pro- 
portional to yt and yb- These are diagrams of order y 2 , y 4 and g 2 y 2 . Because the divergences 
of the theory at zero and at finite temperature are identical, the 3D Lagrangian parameters 
are completely finite and independent of renormalization scale or scheme when expressed in 
terms of the physical observables. 

2 Finite-temperature effective lagrangian 

The procedure for reducing (|3|) to the finite temperature 3-D theory is straightforward. One 
starts with the same 1-loop Feynman diagrams as at zero temperature; the graphs relevant 
for the present problem are shown in figures 1 and 2. However, at finite T, the integrals over 
Po become sums over Matsubara frequencies, po — > l-nnT for bosons and po — *> (2n + l)irT 
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Figure 2: The 1PI diagrams needed for the 4-point function. For explanations see figure 1. 

for fermions, and / d 4 ~ 2e p — > 2ttT J2 n d 3 ~ 2e p. The sum goes over all n ^ for the bosons 
and all n for the fermions to obtain the effective 3-D theory of the zero Matsubara frequency 
modes of the Higgs and gauge bosons, eq. (|l|). Just as for T = 0, one can use dimensional 
regularization (or dimensional reduction in the case of SUSY) to regulate the ultraviolet 
divergences; the divergent counterterms are exactly the same for T > as for T = 0. 
Defining rj tjb = 3^ 2 6 /167r 2 , L b = lnQ 2 /(4vrT) 2 + 2 lE and L F = IuQ 2 /(-kT) 2 + 2-f E , where 
7£ ~ 0.5772 and Q is the arbitrary renormalization scale, the resulting finite-T effective 
Lagrangian is 

= Aree + rj t L F \dH 2 \ 2 + r] b L F \ dH x | 2 
+ (f y 2 T 2 - VtL B (m 2 Q + ml)) \H 2 \ 2 - T} t L B \iiH x + A t H 2 \ 2 



C 3 /T 



+ 
+ 



y 2 T 2 - rj b L B (m 2 Q + m 2 D )) \Hi\ 2 - r] b L B \fxH 2 + A b Hi 



3„,2rp2 
4 
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'' L B (m 2 D + m 2 Q ~2m 2 u )(\H 1 \ 2 -\H 2 \ 2 ) 



32tt 2 

.2 



Vt y 2 (L F - L B )\H 2 \ 4 - \ Vt (g 2 + g'^L^H^ 2 ^ 2 - \H 2 \ 



+ r} b y 2 b (L F - Lb)^ + \ Vb {g 2 + g'^L^H^ 2 ^ 2 - \H 2 \ 



(4) 



where we have used the MS subtraction scheme, which means that the combination 1/e + 
ln47r — 7# has been subtracted from divergences. However the factors of ln47r and je arise 
in a different way at finite temperature, which is why they still appear in the quantities L B 
and L F . Eq. (fj) is an expansion in mg ^/T 2 and we have accordingly dropped all terms 
of order (^/T) 2 and (A t , b /T) 2 . 

This is not yet the complete result for £3; so far we have only integrated out the nonzero 
Matsubara frequency modes, which have masses of the order ttuT. However there still 
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remain particles with masses intermediate between this "superheavy" scale, and the light 
scale which is of order the magnetic mass of the transverse gauge bosons (g 2 T). So we have 
to also integrate out the zero-Matsubara- frequency modes (called "heavy") of the squarks, 
gauge bosons, and Higgs bosons, to an accuracy of y 2 , y 4 or g 2 y 2 in £ 3 . The diagrams that 
contribute are identical to the ones we already considered in deriving eq. (|]); the difference 
is that the heavy particle masses are no longer given by 2imT, but rather the Debye mass. 

To integrate out these remaining heavy modes, we must therefore determine their Debye 
masses, which consist of a tree-level part plus a thermal correction. For the left- and right- 



handed third-generation squarks in a vanishing background field (Hi = 0) one finds |T0 



m 2 = m 2 + ^_ T 2 + yl±y! T 2 + 9^ T 2 + f_ T 2 



'"■ ' '- '9 '6 '4 108" 

m 2 = 2 s + M T 2 + y[ T 2 + V! T 2. 

tR u 9 3 27 

ml = ml + ^T 2 + y lT 2 + 9 ^T 2 , (5) 

where the contributions from gauginos and charginos have been omitted, under the implicit 
assumption that they are so heavy that they decouple. (We have checked that including 
the latter contributions in the Debye masses has no qualitative effect on our subsequent 
numerical results.) In fact these are the only Debye masses we need because the squarks are 
the only particles in L 3 whose tree-level couplings are proportional to y or y 2 . Ignoring the 
thermal loop corrections to these couplings, since they would only give two-loop corrections 
to A3, the result of integrating out the heavy modes of the squarks is 

%\^H 1 + A tdl H 2 \ 2 + ^ Db 

C t M 2 Dt S\ - ( b M 2 Db S b + - DrdH.l 2 - \H 2 \ 2 ) 



C 3 /T = C 3 /T + -^ r \ f id l H 1 + A t d l H 2 \ 2 + -^\fid l H 2 + A b d l H 1 

OlVI rv OlVl r 



Cty 2 M 2 Dt _ ( b y 2 M 2 Db 

4 mr mz ~ 4 mr mr 

( t g 2 M Dt / , Ag 12 f mr 1 



m t L V 3 9 \ m 't R 4 
m b L V 3 9 2 \m~ bR 2 
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where we defined M Dq = m~ qh + m~ qR , ( q = 3y 2 q T/AnM Dq , S t ± = \H 2 \ 2 ± M^^Hi + A t H 2 \ 2 , 
S b ± = \Hi\ 2 ± M D 2 \fiH 2 + A b Hi\ 2 and 

g' 2 T 

Dt = -^~ ( 2m *« ~ m ^ ~ m • ( 7 ) 
The effective lagrangian so obtained is almost in the desired form, but it still depends on 
two Higgs doublets rather than one. At the phase transition, only one linear combination 
of the two is massless ($z), while the orthogonal direction is a heavy field (<&h) which must 
also be integrated out. But since there are no self-couplings of the Higgs fields proportional 
to quark Yukawa couplings, this final step induces no new terms of the order of y A or y 2 g 2 
in £ 3 ; it is just a matter of projecting out the heavy field. Let the angle a describe the 
direction in field space whose eigenvalue in the temperature-dependent mass matrix of the 
two Higgs fields vanishes, at the phase transition temperature T c : 

( Hi \ = ( cos a -sin a \ / \ 

\H 2 ) \ sin a cos a J \ § h J ' ^ ' 

Then the effect of integrating out the heavy field $^ is simply to replace Hi by cos a $/ and 
if 2 by sin a $/. 

Further let m 2 ieS be the entries of the Higgs mass matrix, analogous to the tree- level m 2 's, 
but now corrected by the thermal loop diagrams. If we define the matrices 

™*=(Zl Zl): A,= f, f t »?A; A»=M] "i'V (9) 



ml ml ) 1 \ jj,A t A 2 J ' I jiA b ji 2 



and 



then m 2 ff can be written as 






(10) 



2 2 1 

m cff = m + 



+ (\y 2 t T 2 - V tL B (m 2 Q + m 2 v ) - ( t M 2 Dt ) P t - ( Vt L B + ( t ) A t 
+ (t -> 6; ml^ m 2 D ) \ 



g' 2 

L B {2m\j - m 2 D - m 2 Q ) + D 2 



V32tt 2 

m 2 + i{5Z T ,m 2 } + n T (0), (11) 
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Figure 3: Diagrams contributing to the effective gauge coupling g% in the heavy scale (zero- 
Matsubara-frequency) integration. 



where in the last line we made some definitions that will be useful later. The mixing angle 
is given by 

2 

~ m l,eS / 10 \ 

Sina=- - A — (.12) 



l.eff 



+ m 



3,eff 



,1/2 



at T c , the temperature at which Det(m^ ff ) = 0. The anticommutators in eq. (|11|) arise due 
to wave function renormalization (rescaling Hi so that the kinetic term in £3 is properly 
normalized) . 

At one loop the gauge coupling g% is not renormalized by the Yukawa couplings. We 
checked this by computing the correction to #3 from the correlator Qi&iAiAi. The four rele- 
vant diagrams are shown in figure 3. After rescaling the fields to the canonical normalization, 
the direct contributions to this correlator are found to cancel those induced by wave function 
renormalization. So even after the heavy scale integration we have the tree level relation 
9l = 9 2 T. 

It is now straightforward to extract from eq. @ the temperature-corrected quartic cou- 
pling of the light doublet $/. We find that 
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Figure 4: The tadpole diagrams needed for the 1-loop effective potential in the present 
approximation. The vertices marked by heavy circles (open circles) come from quartic terms 
of order y 2 (g 2 ), where one external Higgs field is replaced by its VEV. 



Eq. ( |i~3|) gives the number that is directly bounded by the lattice results for the condition 
that the phase transition be sufficiently first order, eq. (fj). However, the renormalization 
scale independence of the Lagrangian (|6]), is not yet apparent, and it contains undetermined 
parameters which must be expressed in terms of physical quantities. Once this is done, £3 
becomes manifestly finite and independent of the scale Q. 

3 Relation to the physical parameters 

The next step is to perform the zero-temperature renormalization to the same accuracy as we 
did at finite-temperature in order to express the parameters appearing in eq. ( |T3"D in terms of 
physical quantities, namely particle masses and the vacuum expectation values (VEV's) of 
the two Higgs fields. The VEV's are determined by minimizing the 1-loop effective potential, 
through the equations 



+ 




(13) 



where = sin 2 a — M D 2 (n cos a + A t sin a) 2 and S b a 



cos 2 a — M D 2 (n sin a + Ab cos a) 2 . 



ml + tan (3 m\ + 



m\ + cot P m\ — 




1 dVi-ioop 

1 dyi_i 00 p 
2v 2 dh 2 



(14) 
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Figure 5: Diagrams contributing to the irreducible two point functions of the higgs fields in 
the broken phase. Vertices induced by spontaneous symmetry breaking do not contribute 
to the CP-odd two-point function. Again the indices on the external legs correspond to the 
top sector. 

where we have split the Higgs fields into CP-even and odd parts, Hi = hi + and Vi = 
(Hi) = (hi). The ratio of VEV's is tan/5 = vijv\. Because we have not yet accounted for 
wave function renormalization at one loop, the v^s are not the physical VEV's, defined to 
be Vi = (hi), but the two sets of fields are related by matrix equations 

h = {Z h ) l l 2 h - (1 + l6Z h )h; x = (Z X ) 1/2 X = (1 + \5Z X ) X . (15) 

The matrices 5Z can also be expressed in terms of the derivative of the 1-loop vacuum 
polarizations of the fields, SZ = (dn(p 2 )/dp 2 )| p 2 =0 = IT(0). 

The minimization conditions (|14]) give us two equations for the three parameters m 2 . 
To determine the third we must compute the physical mass of one of the Higgs bosons. A 
convenient choice is the CP-odd scalar, A . Its pole mass, tua, is determined by 

Det(iy xx + n x (m 2 )-m 2 ) = 0, (16) 

where |V XX = d 2 V /dxidxj is the tree-level mass matrix, 

i v 2 . g 2 + g' 2 ,~2 ~2\( i o \ . . 

iYxx = m + - v 3 ) I Q _ 1 I . (17) 

Of course an analogous expression could be used to relate the parameters to the CP-even 
Higgs masses, but the CP-odd sector is simpler because it contains a massless particle, 
the Goldstone boson. Rather than solving for the exact pole mass, we will renormalize 
at p 2 = 0, which means expanding H x (m\) ~ n x (0) + H' x (0)m A , so that H x (m A ) + m 2 A 
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becomes II x (0) — m^Z" 1 . After some manipulations using eq. (|T4| ) to eliminate the second 
term of (0), and using the explicit form of <iVi_i oop /d/ij (see the appendix), eq. ([U]) can be 
written as 

De^-V" 4 Aco^H <18) 

where 

A = -m^ + VtM t F t (Q) + W A b F b {Q), (19) 

m 2 3 is the off-diagonal element of the symmetric matrix m 2 ee m 2 + |{<5Z x ,m 2 }, F t<b (Q) 
are the corrections due to squark loops, 

±m 2 / O 2 \ 

WH ¥<^r( ln i +1 J- (20) 

and the squark masses are given by |IlTJ 



mf ± = \{ml + m 2 Q )+m 2 t + \Ml cos 2/3 



™L = Uml + m 2 Q ) + ml-\M 2 z cos2j3 



± VlK? - ™d - \C b M 2 cos2/3] 2 + m 2 (A + /itan/?) 2 , (21) 



where C t = 1 — | sin 2 6*w and C b — 1 — | sin 2 # W - 



To achieve the simple form (|T8|) for the v4° pole mass condition, we had to introduce the 



shifted angle /3 X , defined by tan/3 x ee v x2 /v x i, where f xi = (Z~ 1 ^ 2 )ijVj. The relation to the 
physical tan/3 is scale- independent: tan/3 x = tan/3 (l — |AZ 1:L + |AZ 22 + AZ 12 cot2/3) , 
where AZ = Z^ — Z x . In the MS subtraction scheme, which we are using throughout, the 
wave function renormalization matrices of the CP-odd and CP-even fields x an d h are 



O 2 

6Z X = -r) t In ^2 Pt - VtH(ml,ml)A t + (t -+ 6); 
= 

m 2 (m 2 +mf_)\ ^ _ mf - m?_ 

— -J 5 P * + 7T m t sm 2 ^ ^ 2 

t-\- t— I t_|_ t_ 

+ Vt sin 2 26> t ( H(m 2 i+ , m?_ ) - ^ ^ ] A f 
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n. ( rri~ + mf m~ — m~ 

+ -^Mf sin 2(3 \ 2 *~ +C t cos 26 t \ 2 1 B, 
24 \ m~ mi mf mj 

\ t_|_ £_ E_j_ t_ 

37774 / wi? — 7nl / ml + m~ 

+ (t->&), (22) 
with P 12 and A t b denned in (J| -[K]) and 



H 2A t ) ' Ab \ ii J ' 

-1 \ ( 2cot/3 -1 

B ' 1 1 2tan/3 J ' B& " 1 -1 



r I 2/jcot/3 A t cot/3-/i\ n = ( ~ 2A b A b tan/5-/i\ , . 
* ~ 1^ Atcotp-fi -2A t )' b ~ { A b tan[3 - fi 2/jtan/? )' 1 ' 

The squark mixing angles are defined by sin 29 t = 2m t (A t + \i cot /3)/(m~ — m| ) and 
008 26*4 = (m 2 / — mg — |C t m| cos 2(3) /{m 2 — m 2 ) (for the bottom squarks let t — > 6, 
m 2 / — > mf, and cos /3 <-> sin /3) and 

__ / 9 9s 1 / m 2 . + m 2 2mim 2 , mi\ . . 

H(m 2 + ,m 2 _) = -— -± ^ - 2 ± " 2 ln-± . (24) 

2 mi — mi \ mi — mi mi-mir mi/ 



Equation fll8|) has a vanishing eigenvalue corresponding to the Goldstone boson, and the 
nonzero eigenvalue is the mass of the A , 

< = (25) 
sm 2(3 X 



Our procedure for computing m\ is somewhat more complicated than that of ref. |TTJ which 



parametrized all the effects of wave function renormalization in scale-dependent VEV's, 
Vi(Q 2 ), whereas we take the VEV's and hence tan/5 to be physical, scale-independent quan- 
tities. Our answer reduces to theirs if we neglect the wave function renormalization effects, 
i.e. by taking (3 X — > (3 and m 2 3 — ► mj. 

The solution fl2"5|) allows the lagrangian mass parameter m\ to be expressed in terms of 
m\. Then, having found m|, the other elements m 2 2 are determined by the minimization 
conditions ([T4]). It is convenient however, to give the results in terms of the scaled quantity 
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rn 



x3 



(see below eq. (|19"D). We find: 

mj(g) = imisin2/? x ( ta ^ 



cot /3 X y 

+ ri t F t (Q)A t + r/ fc F,(Q)A fe + r] t G t (Q)P t + r/ t G b (Q)P 6 



- §Mf(Z>\ cos 2 p - Z h 22 sin 2 /?) (P b Zh ~ PtZ% 2 ) 



+ ri^(g 2 + g ,2 )D(Q)P 3 , (26) 



where 



G q {Q) - E ml (in ^ + l) - 2m 2 (in & + l) (27) 

and the other new function, coming from the tadpoles of the g 2 |g| 2 |if| 2 terms of eq. (|3]), is 
given by 

D{Q) = E < {HQ 2 /ml) + l) - E ™L (MQ7<) + l) 

± ± 

+C t (m 2 Q - m 2 , + \C t M% cos 2/3)F t (Q) 

-C b (m 2 Q -m 2 D - \C b M 2 cos2(3)F b (Q). (28) 

These are all the expressions needed to solve for the entries of the effective mass matrix 
m 2 ff . It is convenient to do so in terms of the scaled matrix m x , using eq. (f|): 

m 2 ff = m 2 + \{bZ T - 5Z X , m 2 } + n T (0), (29) 



where 5Zt and IXr(O) were defined in (|TT]) . The difference between the wave function 
renormalization matrices at finite and at zero temperature is finite and scale-independent, 
as it must be. Moreover, the Q-dependence appearing in m 2 ( fZB"| ) is precisely what is needed 
to cancel that coming from IIt(O) at the accuracy to which we are working, so that the 
matrix m 2 ff and hence the angle cos 2 a are scale-independent. The result is 



m 



■cff 



= m 2 (T) + U2 Vt c B - ( t )A t + {\y 2 T 2 + 2 Vt c B (m 2 Q + m 2 v ) - ( t M 2 Dt ) P t 
+ \ (^K' m '-)-^) {A t ,m 2 } + ^(2c F + ln^ {P,,m 2 } 



g' 2 



+ (t ^ 6 . m ^ m ^ + [j^cb (2mij -m 2 D -m 2 Q )-D T JP 3 , (30) 
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where the matrix m 2 (T) is as in ( p6| ) but with Q = T, and we define cb = ln47r — 7^ and 
cp = hi7r — j E . To one-loop accuracy, it suffices to use the tree level expression for m 2 
appearing in the anticommutators, which in terms of physical parameters is given by 

< ee = M sin 2(3 ( ^ ~^ ) - \M\ cos 2(3 P 3 . (31) 

In the next section we will identify the regions of parameter space where baryogenesis 
is allowed. One must check whether these parameters are in fact compatible with other 
constraints, including the experimental lower limit on the mass of the lightest Higgs boson 
h°, whose tree level expression vanishes when tan/3 = 1. Since we will find that low values 
of tan (3 are relevant for electroweak baryogenesis, it is important to include loop corrections 
to rriho, which can be large. The pole mass of the h° is determined similarly to eq. (|T6|): 

DetdYhfc - H h (m 2 h ) - m 2 h ) = 0, (32) 

where 

, y ^ m2 + ^(35 2 -5l 3 -f4). ,33) 

Expanding the 2-point function IT/j(m 2 ) and performing other manipulations similar to the 
CP-odd case, eq. ([32|) can be put into the form 

Det(M 2 (/?)-m 2 ) = 0, (34) 



where 



+ f sin 2 29 t I(m\ jm\_ ) A t + 2i] t m\ In t+ 4 l ~ P t 

771 ■ 



with 



2 

Tfl~ 

+ 7] t m t sin 29 t In -±P At + (\M 2 Z sin 2(3 5Z% 2 + U\{Q))B t 

+ U\C t + (*-&), (35) 



I(r) = 2- ^Ifiir (36) 
r — 1 
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and 



n?(Q) = -^Mf sin 2/3 In— ^ + -C,cos2^1n-|± 

2 \ ■m- q+ m- q _ 2 m|_ 



^ = J^(^ + ^ 2 ) sin2^fln^-C t cos2^/(m| + /mf_)|. (37) 

It is easy to see that the Q-dependence in ([35]) due to H\(Q) (H^Q)) exactly cancels that 
of SZ22 (SZ^). We verified that the solution of eq. (@), giving m 2 as a function of tan/3, 



m 



2 



agrees with previous results [|T2 

4 Results 

Having completely determined the quantity X^/gl in terms of the physical parameters, we 
now turn to the question of whether the baryons created during the electroweak phase transi- 
tion can be preserved from subsequent sphaleron interactions within the MSSM. Specifically, 
for what values of the physical parameters does Xs/g 2 satisfy the bound (Q)? A priori there 
seems to be no simple, analytic answer to this question. No single term among the many 
loop corrections in eq. (|Tl~D could be identified as being dominant over the others, and fur- 
thermore the parameter space is large: tan/3, m A o, 171%, mfj, m 2 D , /1, A t , A b . We therefore 
chose to do a Monte Carlo search of this space for values which satisfy the constraint @. 

The preceding list of parameters does not include the gauge or Yukawa coupling constants 
because these can be defined through the tree level relations, (g 2 + g l2 )/g 2 = m 2 z /m 2 v , 
y 2 = m 2 /v 2 sm 2 p and y 2 = m\j 'v 2 cos 2 (3. For the gauge couplings this is a consequence of 
the vanishing of the Yukawa contributions to the beta function at one loop, and our neglect 
of the difference between the pole masses and the masses defined at zero external momentum. 
As for the Yukawa couplings, these appear only in loop corrections for us, so to one-loop 
accuracy it is consistent to use their tree-level values. 

For the Monte Carlo search we found it convenient to take as independent parameters 
those listed in Table I, which shows the ranges over which they were varied. The massive 
parameters were allowed to be as large as 1 TeV. These MC-generated sets were subjected 
to various constraints. The lower limits for the pseudoscalar mass m^o and the squark 
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Figure 6: The projected distributions of the data satisfying the sphaleron washout bound 
(0) obtained from the Monte Carlo run described in the text. Units are GeV. 
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tan (3 




2 

m Q,U,D 


p,A t ,A b 


0.4 
10 


20 GeV 
1 TeV 




1 TeV 2 


-1 TeV 
1 TeV 



Table 1: Minimum and maximum values used in the Monte Carlo of the parameters. 



masses were taken to be 20 GeV and 45 GeV, respectively, and we used the recent top 
quark mass measurement of m t = 175 GeV |13[. We further required that the squark 
masses and mixings be consistent with deviations of the p parameter from unity of less 
than 0.01 (see ref. Specifically, the contributions from the top-bottom quark and squark 
splittings as computed in ref. [TJJ , including squark mixing effects, were constrained to satisfy 
Ap(t, b)+Ap(t, b) < 0.01. Finally, the accepted data were required to satisfy the baryogenesis 
constraint (|2]) and stability of the potential (A3 > 0). The distributions of the parameters 
within this set, as well as histograms of derived quantities like the squark and Higgs masses, 
X^loii an d the critical temperature, are shown in figure 6. As a rough indication of how 
special the 5,600 accepted sets were among all possibilities, including some with unphysical 
masses or couplings, we needed 40 million trials to generate them. The baryogenesis-allowed 
cases thus constitute approximately 0.014 percent of the full parameter space. 

One of the most striking features of the distributions is that tan/? is sharply peaked 
near unity and falls to a small but constant value of approximately 10~ 2 of the maximum 
frequency. Thus it would appear that small values of tan/3 are strongly favored by our 
results. This is somewhat misleading however, because there is a strong correlation between 
tan/3 and m^o, as shown in figure 7. The probability of getting large values of tan/3 is very 
much dependent upon m^o. If for whatever reason it became known that m^o was in the 
region of 40 — 120 GeV, the distribution of tan (3 would be much flatter than is shown in 
figure 6, with very large values being almost as likely as small ones. This correlation is also 
evident in the distribution for m^o, which jumps up at small values, due to the enlargement 
of allowed parameters in the direction of increasing tan f3. 

Moreover the distributions for p and A t show a very clear preference for large mixings. 
There are allowed parameters also close to p, A t = 0, but those corresponding to large 



mixing are much more numerous. The large squark mixing angles are also correlated with 
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Figure 7: The baryogenesis-allowed points in the tan/3-m^o plane. The sphaleron bound 
(0) pushes the solutions close to axes; the region away from the axes is populated by points 
with A3/P3 3> 0.04, in violation of the bound. Units of m^o are GeV. 



the rather large average value of 300 GeV of the critical temperature. Due to the smallness of 
the coupling the bottom squark sector corrections are small, which shows in the flatness 
of the ^-distribution. 

The other variables also display certain preferences. The lighter squark masses peak at 
260 GeV and 480 GeV, showing some preference for a moderately light top squark, though 
not as light as that advocated by the perturbative study in reference ||. Furthermore the 
probability for mfj to be negative, suggested in ref. |7| as an optimum choice for strengthening 
the phase transition, appears to be quite small. The mass of the lightest Higgs boson, m^o, 
is in the range of 40 GeV (the experimental lower limit we imposed) to 130 GeV. The Ap- 
distribution |TJJ sharply increases for large Ap, showing the severity of this constraint on the 
data. It is important to note, however, that large squark mixing angles make it much easier 
to satisfy the p-parameter constraint. Finally, the parameter characterizing the strength of 
the phase transition, X^/g^, is monotonically increasing, reflecting the difficulty of getting a 
strongly first order transition. 

There are two further constraints that we imposed on our accepted data sets. First, care 
must be taken to ensure that the finite-temperature perturbative expansion is not breaking 
down where we need it. In particular, the thermal squark loop contributions as written in 
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Figure 8: The density distributions in (mg / tiT c , mu j 'iiT c ) and {tuq / 'ttT c , mu /irT c ) planes. 
Dotted lines show the borders of the accepted regions rnQ jUD /7rT c < 1. 

eq. (f|) are only correct if the zero-Matsubara-frequency ("heavy") modes are much lighter 
than the nonzero ones ("superheavy"), which requires that mQ t u,D < 2nT. We checked that 
the squark contribution to the quartic piece of £3 in eq. (Hj) differs from the exact result 
(not expanded in m 2 /T 2 ) by less than 5% for iuq ud < and by a factor of ~ 2 at 
f^Q,u,D ~ 2ttT. We imposed the more stringent cut of tjiq^u^d < nT c on our data, which 
reduced the size of the sample roughly by a factor of six. The effect of these cuts is clearly 
seen in figure 8, where we show the density plots of the distributions mQ !UjD /irT c , based on 
a separate run with the constraint tuq^d < 27rT c . 

Second, one has to insure the validity of the heavy scale perturbation theory. A typical 
expansion parameter for integrating out the heavy modes is ( t — ty^T '/ '^kMjj^ which should 
not become too large; we required that T/m < 1 for all Debye masses m. From the distribu- 
tion for T/rrif in figure 6 one sees however that the other cuts (in particular the requirement 
m Q,u,D < tt^c) already confine the heavy scale expansion parameters to the range needed 
to insure ( < 1. This is partly because we did not consider negative values for mfj in the 
present work. 

Clearly, making these consistency cuts can lead us to neglect parameter values that might 
actually be acceptable for electroweak baryogenesis. However it is difficult to consistently 
implement the dimensional reduction program except for the relatively light squark masses 
that pass our cuts, or in the opposite limit of mq t u,D 3> 27rT c where the heavy modes 
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decouple, since only then is there a clear hierarchy between the superheavy and heavy scales. 
(Of course the same problem exists in the purely perturbative effective potential approach.) 
A naive way to interpolate between these limits would be to replace the m 2 /T 2 expansions 
with the exact expressions for the corresponding finite temperature integrals. While not quite 
rigorous, this might provide a reasonable approximation to the exact result. An investigation 
along these lines is in progress. 

5 Conclusions 

We examined the strength of the first order electroweak phase transition in the MSSM with 
respect to the prospects for safeguarding electroweak baryogenesis from washout by residual 
sphaleron interactions, finding rather encouraging results. Although our calculations were 
perturbative, the method — computing the three-dimensional finite-temperature effective ac- 
tion of the light Higgs and gauge fields at the phase transition — enabled us to take advantage 
of nonperturbative results that have been obtained from lattice gauge theory computations. 
This is therefore the first study of the phase transition in the MSSM that can claim to be free 
from the infrared divergences that make the usual perturbative calculations untrustworthy. 
Indeed, we find that the results of these two distinct approaches disagree. 

The most serious constraint on electroweak baryogenesis in the MSSM, seen also in the 
earlier investigations ||-[f§, seems to be the biasQ toward small values of tan/3 < 2. But in 
our approach this bias is very different from a prohibition on large values of tan /3; in fact 
such values are not ruled out, but they must appear in conjunction with low values (40 — 120 
GeV) of the A boson mass, as is clear from our figure (7). Since there is no intrinsically 
"correct" integration measure for the space of all parameters in the MSSM, the large tan f3 
possibility can hardly be considered less natural, even if it comprised a smaller subset of our 
Monte Carlo results. If the limit on m^o should be improved such as to exclude this region, 
then it will become a question of whether tan/3 <, 1.5 is viable. 

There is one important caveat to our conclusions: as emphasized in || , the approach used 
here assumes that the only light degrees of freedom at the phase transition are the transverse 

1 This restriction is somewhat alleviated in a two-loop computation, however [fl5|. 
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gauge bosons and a single linear combination of the two Higgs fields. It is possible that other 
fields which are generically heavy happen to also be light for special parameter values-for 
example, the Debye masses of the squarks can vanish if tuq ud < 0. In such cases we can 
say nothing until the lattice computations are redone to take into account these potential 
new sources of infrared divergences. 

Note Added 

Soon after this paper was completed, we received two articles where the same problem was 
considered. Losada |l6j derived the dimensionally reduced lagrangian, keeping also the g 4 - 



corrections (but setting g' = 0). Laine |17| made a complete analysis in an approximation 



similar to ours. Where comparison is possible, our results are in good agreement. 
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A Appendix 

Here we give the derivatives of the effective potential and the zero-momentum limits of 
2-point functions, which were necessary to derive, but not to finally express, the relevant 
quantities in the body of the text. 
1 dVx 



2v% dhi 

1 dVi 

2v<? dh? 



-7] t (/j 2 + fiA t taa P)F t (Q) - Vb(A 2 b + ^A b tan P)F b (Q) 
-V b G b (Q)-ji^(g 2 + g' 2 )D(Q) 
-r] t (A 2 + \l A t cot /3)F t (Q) - r/ t (/i 2 + fiA b cot (3) F b (Q) 
-VtG t (Q) + ji^(g 2 + 9' 2 )D(Q), (38) 

with F q (Q) defined in (^), G q (Q) in (^7\) and D(Q) in The two-point function of the 
CP-odd sector at the zero momentum is given by 

11(0)* = - Vt G t (Q)P t - r] t F t {Q)A t + (t ^ b) 



20 



128 



^(g 2 + g' 2 )D(Q)P 3 , (39) 



with Pj denned in (|Tl]). In the CP-even sector we find the result 

n(o) h = n(o) x 



2 

+ ^ sin 2 26 t I {ml /ml ) A t + r) t m t sin 26 t In ^P At 

2 + m\ 

2 2 
TTZ- 777.- 

+ 2 Vt ml In f+ 4 - P t + n*(Q)B t + II* C t 

+ (*-&), (40) 



with Pa 9 , Bj and C, defined in (p3|), /(r) in (|36| ) and Ilf in 
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